GLOBAL EXISTENCE OF NULL-FORM WAVE 
EQUATIONS ON SMALL ASYPTOTICALLY EUCLIDEAN 

MANIFOLDS 

CHENGBO WANG AND XIN YU 

^ ' Abstract. We prove the global existence of the small solutions to the 

I Cauchy problem for quasilinear wave equations satisfying the null con- 

dition on (R'',0), where the metric g is a small perturbation of the flat 
metric and approaches the Euclidean metric like (1 + \x\^)~''^'^ with 
p > 1. Global and almost global existence for systems without the null 
condition are also discussed for certain small time-dependent perturba- 
tions of the flat metric in the appendix. 
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1. Introduction 



This paper is concerned with the Cauchy problem for the system of the 
quasi-hnear wave equations in three space dimensions of the form 



subject to suitably small initial conditions, posed on certain asymptotically 
Euclidean manifolds (M^,g). We assume that the propagation speeds, cj > 
0, are distinct, and we refer to this situation as the nonrelativistic case. 

We shall construct a unique global classical solution, provided that the 
coefficients of the nonlinear terms satisfy the null condition and the metric 
g is a small perturbation of the flat Euclidean metric and approaches the 
Euclidean metric like (1 + with p > 1. 

In the Minkowski space-time (with flat metric g), this problem has been 
extensively studied. In the relativistic case (with same c/), the null con- 
dition was first identified by Klainerman and shown to have global exis- 
tence of small solutions in Christodoulou [6] and Klainerman [29] (see also 
Hormander [11] and John [16]). Without the null condition, we have almost 
global existence [18], [28], [16], but as examples show, arbitrarily small initial 
conditions can blow up in finite time for some equations [15], [44]. Notice 
that this does not mean that the null condition is the necessary condition 
for the general quadratic quasilinear problems to have global existence with 
small data. Actually, there is a larger class of nonlinearities which will en- 
sure global existence, which is related to the so-called "weak null condition" 
of Lindblad and Rodnianski, see e.g. [35], [36], [4], [24]. 

Small solutions always have global existence in higher dimensions [31], 
[43], [28]. The two-dimensional relativistic case is rather more complicated. 
The sharp results are given in [2], [3], with previous works in [12], [19]. 

In the non-relativistic case, the null condition guarantees that the self- 
interaction of each wave family is nonresonant and is the natural one for sys- 
tems of quasilinear wave equations with multiple speeds. It is equivalent to 
the requirement that no plane wave solution of the system is genuinely non- 
linear, which follows from the so-called "John-Shatah's observation" (John 
[17] (p. 23) and [1]). Related results were established, e.g., by Kovalyov 
[34], Agemi and Yokoyama [1], Hoshiga and Kubo [14], Yokoyama [52], 
Sideris and Tu [46], Kubota and Yokoyama [33], Sogge [47], [48], Hidano 
[8], Katayama [20] [21], Katayama and Yokoyama [25], and Katayama and 
Kubo [22]. It is remarkable that the approach of [46] and [8] does not use 
estimates of the fundamental solution for the free wave equation, which is 
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more robust when considering problems with variable coefficients. The two- 
dimensional case has also been examined, see [13], [14], [8] and references 
therein. 

In exterior domains, null form quasilinear wave equations were previously 
studied by Keel, Smith and Sogge [26], Metcalfe and Sogge [40] [42], Met- 
calfe, Nakamura and Sogge [38] [39], and Katayama and Kubo [23]. The 
general quasilinear problems were also studied, see [27], [41] and references 
therein. 

It is interesting to investigate similar problems on various space-time man- 
ifolds. Recently, there have been some progress in this direction. Global and 
almost global existence of the semilinear problems posed on asymptotically 
Euclidean manifolds have been obtained in Bony and Hafner [5] , Sogge and 
the first author [49] , and the authors [50] . Global existence for the null form 
semilinear wave equations on slowly rotating Kerr spacetimes or time de- 
pendent inhomogeneous media has been obtained in Luk [37] and Yang [51]. 
In this paper, we will deal with the nonrelativistic quasilinear wave equa- 
tions on small asymptotically Euclidean manifolds, mainly inspired by the 
approach of [46] and [8] , together with the recent local energy estimates with 
variable coefficients obtained in Hidano, the first author and Yokoyama [9] . 
It will also be very interesting if we can deal with the problems with general 
non-trapping asymptotically Euclidean manifolds or small asymptotically 
flat manifolds with certain time-dependent metrics. 

Before stating our main result, we introduce the necessary notations. 
Points in will be denoted by {x^ ) = (t,x). Partial derivatives 

will be written as da = d/dx°', with the abbreviations d = {80,81,82,03) = 
(9t,V). Here, we have used the convention that Greek indices range from 
to 3 and Latin indices from 1 to 3. We will also abuse the notation to 
use the Greek indices to denote multi- indices, which should be clear in the 
context. Hereafter, the Einstein summation convention will be performed 
over repeated indices. The angular-momentum operators are defined as 

n = (01,172,03) = xAV, 

where A denotes the usual vector cross product in M^, and the scaling oper- 
ator is defined by 

(1.2) S = t8t + r8r = x'^da, r = \x\. 

The collection of these eight vector fields will be labeled as 

r = (ro,...,rr) = {8,n,s). 

We consider the asymptotically Euclidean Riemannian manifolds (M^,0) 
with 

= gij (x) dx^ dx^ . 

The metric g is assumed to be a small perturbation of the flat metric. More 
precisely, we suppose gij{x) E C°°(M^) and, for some fixed p > and 5 <^1, 

(HI) Va G \d^{gij - 5ij)\ < 
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with 6ij = 6^^ being the Kronecker delta function and (x) = -^/i + |xp. 
Since (5 ^ 1, it is clear that the metric g is a non-trapping perturbation. Let 
g = det{gij), the Laplace-Beltrami operator associated with q is given by 

where {g^^ (x)) denotes the inverse matrix of {gij{x)). 

Consider the initial value problem for the nonlinear equations of the form 

(1.3) {U^uY = {d^t -c]\)u' = N\u,u),I = 1,2,- ■■ ,M 
in which the quadratic nonlinearity = Q + 5" is of the form 

(1.4) Q\u, v) = Qjf^d^u-^d/.d^v'', S'iu, v) = S'jfd^u-^dpv''. 

The construction of solutions will depend on the energy integral method, 
which requires the quasilinear part to be symmetric: 

(1.5) qy/-^ = qy^' = Qfr^'. 

The key assumption for global existence is the following null condition which 
says that the self- interaction of each wave family is nonresonant: 

(1.6) Q'jf^am^ = S\f^^ip = for all ^ s.t. = + ^2 + Cl) • 
The standard energy norm is denoted as 

Ei{u{t)) = -Y, / \du\t,x)\''dx, 

and higher order derivatives will be estimated through 

(1.7) EMt))= EiiT'^uit)), m = 2,3,--- . 

a|<m— 1 

In order to describe the solution space, we introduce the time-independent 
vector fields A = (Ai, . . . , A7) = (V, il, rdr). Define 

H^{R^) = {/ G l2(IR3) . A"/ G L^, \a\ < m}, 

with the norm 

(1-8) \\f\\H- = E ||A"/||l2. 

|ci| <m 

Solutions will be constructed in the space H^{T) obtained by closing the 
set C7°°([0,T);C^(M3)) in the norm sup EU^{u{t)). Thus, 

0<t<T 

By Sobolev embedding, it fohows that Hp{T) C C™-2([0, T) x 
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An important intermediate role will be played by the following two weighted 



norms 



(1.9) Xmiu{t)) 



M 

EE E IKc/t-r)5-rVwiiL^ 

1=1 \a\=2 \/3\<m-2 



and 
(1.10) 



3m{u{t)) 



M 

E E 

1=1 |a|<m-l 



L2 



with /X € (0,1/2) and fi' > fj, to be determined later (the choice will be 
II = 1/4 and /i' = min(2p — l,3)/4, see (3.1), (3.5)). The second norm is 
extracted from the local energy norm (also known as KSS-type estimates), 
which is defined as 

(1.11) LErn{t) = [ 3rrMr))dT . 

Jo 

Let us now state our main result precisely. 

Theorem 1.1. Let p > 1 and 6 <^ 1. Assume that the nonlinear terms 
in (1.4) satisfy the symmetric and null conditions (1.5), (1.6). Then there 
exist constants e ^ 1 ^ Co, such that the Cauchy problem for (1.3) has a 
unique global solution u £ H^{T) for every T > 0, when the initial data 

du{0) G H1~^{R^), K>9 

satisfying 

(1.12) E'J^,iu{0)) exp CoEl/\um < e. 
Moreover, the solution satisfies the bounds for some Ci, C2 > 1, 
LE^^2{t)+E^^2{u{t)) < 2Cie and LE^{t)+E^{u{t)) < C2E^{u{0)){tf^' . 

To conclude the introduction, let us give some remarks and comments. 

Remark 1. We remark here that, the situation for the null form quasilinear 
problems seems much more delicate than the general quasilinear problem, 
technically due to the occurrence of the scaling vector field in our argument. 
Actually, it is not hard to see that we can prove the almost global exis- 
tence (and global existence for higher dimension) for the solutions to the 
quasilinear quadratic equations, on asymptotically flat manifolds with small 
time-dependent metric perturbation, by combining our argument with the 
approach of Metcalfe and Sogge [41] , where there is no need to use the scal- 
ing vector field. See Appendix 5 for the proof. Although there is work [22] 
dealing with null form problems without using the scaling vector field in the 
literature, it seems difficult to be adapted for the setting of time-dependent 
asymptotically flat perturbation (except the case of compact perturbation, 
see [23]). 
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Remark 2. The case n = 2 seems more difficult to handle, mainly because of 
the lack of the local energy estimates with variable coefficients. In principle, 
we should be able to prove global existence under the null condition of 
Alinhac [2], provided that we have such local energy estimates, by combining 
our argument and that of Hidano [8] . 

Remark 3. The same argument can yield global results for the system with 
repeated speeds (including the relativistic case), by strengthening the null 
condition to be nonresonant interaction among the waves with the same 
wave speeds. See [46] or Section 5 of Chapter II of [48]. 

Remark 4. One may wonder if the standard local energy estimates with 
variable coefficients in [41] (with = 1/2 instead of 1/4) is sufficient in our 
argument. It seems to us that it is not the case, see e.g. (3.2) and (3.5). In 
this sense, it is the work of [9] which makes the adaption of Sideris and Tu's 
argument possible. 

2. Preliminaries 

2.1. Commutation and null forms. In preparation for the energy esti- 
mates, we need to consider the commutation properties of the vector fields 
r with respect to the nonlinear terms. It is necessary to verify that the null 
structure is preserved upon differentiation, in some sense. 

Lemma 2.1. Let u be a solution of (1.3). Assume that the null condition 
(1.6) holds for the nonlinearity in (1.4). Then for any a G N®, 

^+7+At=o |/9|<|a|-l 

in which each is a quadratic nonlinearity of the form (1.4) satisfying 
(1.6), and with m £ N denote functions such that 

|V"r^| < C„5(r)-''-™"l"l for any a G . 

Moreover, if = 0, then N^^ = N . 

Proof. Define 

[r, N]{u, v) = TN{u, v) - N(Tu, v) - iV(n, Tv). 

This is a quadratic nonlinearity of the form (1.4). Moreover, by the proof 
of Lemma 4.1 in [46], we know that [T, N] is null form for each T. 
By (HI), we have 

□g = □ - roV^ - nV , 
with □ = Diag{dj - A, • • • , - 4^A). Recall that 

we want to prove the result by induction. It is clear the result is true for 
|a| =0. Now assume that it is true for any a with |a| = m. Given ag with 
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laol = m + 1, we can find some j and a with \a\ = m and F"" = r^T". 
Then by the inductional assumption, we can calculate as follows 

□gr"Ou = U^TjT'^u 

= 25j'jUT'^u - [roV^ + n V, Tjlr^u 

^+7+M=« |/3|<|«|-1 
l/3|<l"l 

+ {[r,,A^;](rVr^n) + iv;(r,rVr^«) 
+iv;(r^n,r,nn)} 

which is of the required form. This completes the proof. □ 



2.2. Estimates for null forms. The utility of the null condition is cap- 
tured in the next lemma, where we get some additional decay in nonlinear- 
ities with the null structure (1.6). 

Lemma 2.2. Suppose that the nonlinear form N(u,v) defined (1.4) satisfies 
the null condition (1-6). For any u, v, w £ C^([0, T] x M^) and r > C/t/2, 
we have at any point (t, x) G [0, T] x M'^ 



(2.1a) \Q 



C 



II 



< 



(r)L 



Tu\\d v\ + \du\\drv\ + \du\\dv\ + {cjt — r)\du\\d v\ 



and 



^/,o/37^ 

< 



(2.1b) \Q{i daudfsvd^w 
C_ 

V) 



Tu\\dv\\dw\ + \du\\Tv\\dw\ + \du\\dv\\rw\ 

+ {cjt — r)\du\\dv\\dw\ 



(2.1c) \Sjf'^daudpv\ < -^yru||9i;| + \du\\rv\ + {cit - r)\du\\dv\\dw\ 



C r 



Proof. The inequalities (2.1a)-(2.1b) are exactly Lemma 5.1 of [46]. The 
proof of (2.1c) is similar. See also Lemma 5.4 of Chapter II in [48]. □ 
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2.3. Sobolev-type inequalities. The following Sobolev inequalities do not 
involve the Lorentz boost operators. The weight (c/t — r) compensates for 
this. We use the notation defined in (1.7), (1.9). 

Lemma 2.3. We have the following inequalities for smooth functions u : 

Tn>3+1 , -mM 



(2.2) 
(2.3) 
(2.4) 
(2.5) 
(2.6) 



^/^\u{t,x)\ < CEl^^uit)), 
\du{t,x)\ < CEl^^{u{t)), 

{cjt-r)^/^\du^{t,x)\ < c\El/^{u{t))+X3{u{t)) 
{cit - r)|aV(t, x)\ < CXi{u{t)), 
^^^{cit-r)\du^{t,x)\ < c\El/\u{t)) + Xs{u{t)) 



See Lemma 6.1 in [46] for the proof of (2.2)-(2.5) (see also [30], [45]). The 
inequality (2.6) is just (4.22) of [8]. See also Lemma 5.2 of Chapter II in 
[48]. 

Let us give the proof of (2.6) here. Recall the following trace estimate 

r-'/'ll/(™)||L4 <C||/||^, 

(see (3.16) of Sideris [45] and also (1.3) of [7] for more general version), then 

r^/^{cit-r)\du^{t,x)\ < C ^ r^/^\\{cit-r)n"du^{t,ruj)\\L4 

|q|<1 

< \\{cit-r)n''du\t, 



|q|<1 

< c[El'\u{t)) + X^{u{t)) , 

which is exactly (2.6) for r > 1. For the remained case r < 1, we can use 
a cut-off argument to derive the required estimate as follows. Let x be a 
cut-off function such that X = 1 for |x| < 1 and x(^) = for |x| > 2, then 
for r < 1, 

(r)i/2||(c/t -r)9n^ (t,x)||ioc(|,|<i) 

< C{t)\\xdu\t,x)\\L^ 

< C{t)\\xdu'{t,x)\\j,2 

< C{t) (^5^||vW(t,x)||i2(|,|<2) + ||9n^(t 

2 

< CY,\\{cit - r)V'du'{t,x)\\L2 + C\\r^/^{cit - r)du\t,x)\\L^ 



4 = 1 



< c 



[El/\u{t))+Xs{u{t)) 
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where we have used (2.6) for r > 1 in the last inequahty. 

2.4. Local energy estimates. One of the main extra steps in our proof 
is to exploit the local energy estimates (also known as Morawetz estimates, 
KSS estimates), to handle the extra terms arising from the non-flat metric. 

Lemma 2.4. Let fo = (r/(l+r))^'', = r/{r+2'^) with k > I, and u he the 

solution to the equation (df - cjA + h^^°'l^{t,x)dadp)u^ = in [0,T] x M" 

with /i^'"^ = /i^'/^", E/liEo<a,/3<nl^^'"^l < 1/2 and n > 3, then there 
exists a constant C > 0, depending only on the dimension n, such that 

sup Ei{u{t)) + LEi{T) + (log(2 + T)Y^KSSi{T) 

0<t<T 



< CEi{u{0)) + C [ [ \dh\\du\'^dxdt 
Jo Jr" 

F^dtu^dxdt 



+C 



M 

E 

1=1 







+C sup 

k>0 



M 

E 

1=1 
+Csup / 

A:>0 JO 

< CEi{u{Q)) + 



fk dru' + 



n 



2r 



-u 



F^dxdt 



fk\dh\ + \fl\ + 



fk 



\h\ 



\du\ \du\ + 



+C 



M 

E 

1=1 



F^dtu^dxdt 



\du\ \du\ + — 



dxdt 
dxdt 



+C sup 

fc>0 





M 

E 

7=1 



T 



fk drU + 



n 



2r 



-u 



F^dxdt 



See Section 2 of Hidano, Wang and Yokoyama [9] for the proof (see also 
[10], [41] and references therein). 

2.5. Weighted decay estimates. One of the main steps is to control the 
weighted norm Xfi{u{t)). This will be accomplished in this subsection by a 
type of bootstrap argument, similar to what in [46]. 

Lemma 2.5 (Klainerman-Sideris estimate). Let u G H^iT), 5 <^ 1 and 
p > 1. Then 



(2.7) 



X^iu{t)) < C 



El'\u{t))+ \\{t + r)n,T^uy-c)u^ 

\a\<K~2 



Proof. The same estimate is known to be true for the standard D'Alembertian 
□ instead of Dg, see Lemma 7.1 in [46], see also [32] and Lemma 5.3 of 
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Chapter II in [48] . To complete the proof, we need only to control the norm 
involving □ by that of Dg and good terms, as follows. For any a E with 

\a\ < K — 2, 

II (t + r) aT^u{t)\\L. -\\{t + t) n,V^u{t)\\L2 

< \\{t + r) r^V^T^u{t)\\L2 + \\{t + r) nVT^u{t)\\L2 

< E (11^^^* - ^) ^-iV'rV(i)||L^ + \\{cit - r) roVrV(t)||^0 

< C5X^{u{t)) + C5\\V{{cit - r)VrV(t))||i2 

< C5X^{u{t)) + C5\\Vn{t)\\L2 

< C6X^{u{t)) + C5El/\u{t)) , 
where we have used the elementary inequality 

(2.8) t + r <C{cit-r){r) 

in the second inequality, the Hardy's inequality and p > 1 in the third 
inequality, and the fact that \V{cit — r)\ < C in the fourth inequality. □ 

Now we assume that u solves the nonlinear equation (1.3). 

Lemma 2.6. Let p > I and u G H^{T) be a solution of (1.3). Define 
k' = [^] + 3. Then for all |q| < k - 2, 

(2.9) \\{t + r)U,T^u{t)\\L2 < C[XMt))E]J\{u{t)) + XMt))E]!,\um 

+ CEf{u{t))E'J^Mt)) + C5X^~i{u{t)) + C5E'J^,{u{t)). 

Proof. There are similar estimates for the case of □ in Lemma 7.2 of [46] 
and Lemma 5.2 of [8]. In view of Lemma 2.1, we need to control the terms 
of the form 

||(t + r)9rV(t)a2rV(t)||i2 , 
||(t + r)arV(t)5rTu-^(t)||i2 , 

with |/3| + I7I < |a| < K - 2 and 

||(t + r)roV'T^u{t)\\l, + ||(i + r)r,Vr^u{t)U. 

with |/3| < |a| — 1 < K — 3. For the first set of terms, we separate two cases: 
either |/3| < — 3 or I7I < k' — 4. In the case of |/3| < k' — 3, using (2.8) 
and (2.3), we have 

\\{t + r)dT>'u'{t)d^r''u-'{t)U2 < C\\{r)dTf'u'{t)\\L^\\{cjt-r)d^r^u'{t)U2 

< CE'f{u{t))X^{u{t)) . 
In the second case, using (2.8) and (2.5), we get 

\\{t + r)dT^u\t)d^ru\t)\\L2 < C\\dr^u\t)\\L2\\{r){cjt-r)d^r"'u-'it)U^ 

< CE'J^Mm^'Ht)) . 
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Turning to the second set of terms, without loss of generaUty, we can 
assume |/3| < I7I (and so < k' — 3). Then, using (2.6) and (2.5), we get 

\\{t + r)dT'^u^{t)dT^u'^{t)\\L2 

< CiWidt - r)arV(t)||L.o + \\{r)dT^u'it)\\Loo)\\dT^u'it)h2 

< C{E'f^Mi)) + X.'HmElL'Mt)) + CE'f{u{t))E'J_\{u{t)) . 
Since /O > 1 and < |a| — 1 = k — 3, we see that 

II (t + r)roV^TPu\t)\\L2 + \\{t + r)riVrV(t)||i2 

< C7|| {at - r)r^^V^TPu\t)U2 + \\ {at - r)roVTPu\t)\\l, 

< C6X^^,{u{t)) + C6\\V{{cit - r)Vr^n(t))||i2 

< C6X^^iiu{t)) + C6El^^^{u{t)) . 

This completes the proof. □ 

The next result gains control of the weighted norm X by the energy. We 
distinguish two different energies, the smaller of which will remain small. In 
Section 4, we will allow the larger energy to grow in time. 

Lemma 2.7. Letu G H^iT), k>8, be a solution of (1.3). Define t] = k—2, 
and assume that 

< 1, eo = sup El/'^{u{t)) < 1, p > 1. 

0<i<T 

Then forO<t <T, 

(2.10a) X^{u{t)) < C^y2(u(i)) 

and 

(2.10b) X^{u{t)) < CEll\u{t)). 

Proof. Let r]' = + 2>, rj = k — 2. Since 77 > 6, we have 77' < t]. Thus, 

by Lemmas 2.5 and 2.6, we get from our assumption 

X,(n(t)) < C[El/\u{t)) + Ell^Xr, + X^,El':^^ + EH''El'\] 

+C5[X.,^^{u{t)) + El%{um 

< C[Ey\u{t)) + {eo + 5)X^{um 

Thus, if eo and 5 are small enough, the bound (2.10a) results. 

Again since k > 8, we have kI = [^^^] + 3<r/ = ft; — 2. From Lemmas 
2.5 and 2.6, we have 

X.{u{t)) < C[Ell\u{t)) + X,E]!:',+X^El'^ + Ell^E'J\] 

+C6[X^^^{u{t))+E]/^^{um 

< C(l + + E]'^ + 5)El'^{u{t)) + C{Ey^ + 5)X^{u{t)) 
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If we apply (2.10a), eg <C 1 and 5 <C 1, then 

^n{u{t)) < CEl/\u{t)) + C(eo + 5)X^{u{t)), 
from which (2.10b) follows. □ 

3. Energy and local energy estimates 

Assume that u{t) S H^{T) is a local solution of the initial value problem 
for (1.3) (in which we need the symmetry condition (1.5)). Our task will be 
to show that E^{u(t)) remains finite for all t > 0. To do so, we will derive 
a pair of coupled integral inequalities for Ei^{u{t)) + LEi^{t) and Ej^{u{t)) + 
LE^(t), with r] = K — 2. If (1.12) holds with 2CiE < Eq occurred in Lemma 

1 /2 1 /2 

2.7, then E"^ (^^(0)) < e and Er^ {u{t)) < 2Ci£ for certain small interval 
t G [0,r]. Define 

To = sup{T : El/^{u{t)) < 2Cie,t G [0,r]} . 

Here the constants Co, Ci > 1 will be determined later. All of the following 
computations will be valid on the time interval [0,To). 

To complete the proof of Theorem 1.1, we need to obtain estimates for 
Em{u{t)), as well as Jmiu{t)) with m = r], k. Instead of giving the estimates 
for 3m{u{t)) directly, we will give the estimates of LEm{u{t)). 

Since we have 

{U,Y = - cjA - c]{g'^ - 6'^)did, + nV , 

it is easy to see that Lemma 2.4 holds for (D)^ — h^'^^didj with /i^'*-' = 
Cjlg"^^ — 6^^). By Lemma 2.1, we see that 



(□gT^)^ - riVr^n^ 



= ^^M"(rVr^n) + 

= Fa + Ga ■ 

Since p > 1, if we assume 



l/3|<l«hl 



(3.1) 



2/i' < /) - 1 + 2/i , 
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then Lemma 2.4 tells us that 
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(3.2) sup Em{u{t)) + LE^{T) 

tg[0,T] 

< CE^{u{Q)) 



\a\<m—l 
|a| <m— 1 



M T 
^ 'O J9? 



1=1 



\a\<m—l 



sup 

fc>0 
rT 



M 



E 

7=1 



fk ( drT'^u' + -Vu' ] F^dxdt 



C Y J , (^\dr''u\ + ^\T''u\^\Ga\dxdt 



\a\<m—l 

< CEM0)) + C5LEm{T) 



+c E 

|a| <7n— 1 



JR^ 



Y F^dtT^u^dxdt 



1=1 



-C 2, sup 

|a|<m— 1 — 



/" /" y /fc f 9r-r V + -r V Vidxdt 

CS;„(n(0)) + C75LE^(r) + C V r CUt)dt 

+C V sup rc^^kit)dt 
1 fc>o Jo 



Q:|<m— 1 



|a|<m— 1 ' 



dxdt 



3.1. Estimate for Cf . Notice that 



where A'^^'" denotes the /-th component of N'j^. Among all of the terms, the 
cases |q;| = I7I = m — 1 for QJJ are quasilinear terms and we want to use the 
symmetry condition (1.5) to absorb such terms. When |a| = I7I = m — 1, 
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we have Q' 
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= Q, = ot. Then an integration by parts argument will yield 



V / Q\u,T'^u)dtT'^u\t)dx 

M „ 



7=1 

M 



1=1 



-d^dpu-^ d^V'u^ dtV'u^ - d^u^df.V'u^d^dtV'^u^] dx 



M 



1=1 



J2 QjE^'dt / dpu^d^r''u''dtT''u'dx 



M 



JK 



1=1 



dx 



1=1 



M 



dx 



where we have used the symmetry of the equation (1.5) in the third step. 
By introducing the notation 



M 

(3.3) Efit) = JZiQjfS^o - i^QjK^) / d^u'd^n^'d^r^u'dx 
1=1 -^^^ 
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we see that 



(3.4) rcfm 



a|<m— 1 
M 

S E E 

1=1 |q| <m— l,/3+7+^=o 
M 

+E E 

1=1 |o[<m— l,/3+7+^=Q,|7[<m— 1 
M „x 

^1 

M „x 



[ [ Sl'"(T'^u,T^u)dtr"u^{t)dxdt 
Jo 

I I Ql:'^{Tf^u,T^u)dtr''u\t)dxdt 
Jo Jrs 

fl Yl 11 QjK"'d-ydi3U-^d^r''u^dtT''u^dxdt 

7=1 -^0 Jr^ 

+ V V / / QYK^dpdtu-^d^T''u^d^T''u^dxdt 

frti I 1 -'0 Jm? 

1—1 \a\=m—l 



+ \EnT)-E'^m. 

|«|=m— 1 



3.2. Estimate for C2 Let = fk{dr+\) = h\{x)di+hj^{x) with /ifc(x) 
and h\{x) = x^hi^{x), we see that 



» M M . 

-/=! 1=1 p+7+^=« 



As for Cf , for the case \a\ = I7I = m — 1, we have Q^°^ = Q, ^ = a, and 

M „ 

y / Q^(u,r°n)Lfcr"n^(t)(i2; 
A/ „ 

= y / QY/''dpu^d^i^''u'^LkT''u'dx 
1=1 J^^ 

M „ 

1=1 -^^^ 

M 
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The last term in the above identity can be estimated as fohows, 



M 



1=1 

M 



1=1 

M 



1=1 -'^^ 

QYk^ I dfsu-^d^n^ [{d-yhi)di + {d^hk) + h^d^] T'^u^dx 



+l:Q 



M 



[{dihi)dpu-^ + hid/sdiu-^] d^V'u^d^r^u^dx 



1=1 

M 



{d^hi)di + [hk- -d,hi ) + d^hk 



V'u^dx 



1=1 ''^^ 



where we have used the symmetry of the equation 1.5. It is easy to check 
that, by setting, say, /i = 1/4, there is a uniform constant C > which is 
independent of > such that 

|V/ii| + \hk\ < Cr-i/2(r)-i/2^ |v/ifc| < C7r-3/2(r)-V2. 
Moreover, by (2.3), we have 

\du\ < C{r)-^El^^{u{t)) < CCie{r)-^ . 
Then, the first term in the above expression is controlled by CCieJrn ^t) with 



(3.5) 



fj, = l/4< n' < 3/4 . 



In summary, by introducing the notation 



M 



(3.6) 



1=1 



I,l3ii0 
JK 



dfiU-^d^V'u^LkT'^u^dx, 
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(3.7) 



we have proven the following 

\a\<m-l •'^ 

M 

S E E 

1=1 |Q|<m-l,/3+7+M= 



E 

|o| <m— l,/3+7+^=o,|7| <-m— 1 



+ 



+ 



/ / Sl;''{T^u,T'^u)Lkn\t)dxdt 
Jo Jrs 

[ [ Q^''^{T^u,T"'u)Lkn\t)dxdt 
Jo Jr^ 

M T r 

Jl 12 Qji'^^'d^disu-^d^n^Lkn^dxdt 
I I , Jo Jr^ 

AI „T „ 

QYK^^kdpdiU-^d^V^u^d^V'u^dxdt 

J. I I , Jo Jr^ 

-1=1 |o|=m— 1 



1=1 |o|='m— 1 

+ E \EUT)-El^m+CC^eLE^{T) . 

\a\=m—l 

4. Proof of the energy inequality 

In this section, we will complete the proof of Theorem 1.1. 
At first, we claim that by (3.2), (3.4) and (3.7), we can deduce that for 
some universal constant C3 > 1, 

(4.1) LE^{T)+ sup E^{u{t)) 



0<t<T 



H0)) + C3 [ {t)-^Ell\u{t))EMt))dt 
Jo 

-CsC,i6 + e)LE^iT) + Cs ^ sup \E^^,iT) - E^^.m 

I 1 k>0 



|q!|=K— 1 



and 
(4.2) 



+C3 E \EnT)-Efm 

\a\=K-l 

LErj{T)+ sup E^^{u{t)) 



0<t<T 



< CsE 



V,(n(0)) + C3 r {trh~^I^El'\u{t))E.,{u{t))dt 
Jo 

+C3Ci{6 + e)LEr,{T) + C3 V snp\E^ ,{T) - E!^ 

I , k>0 



a|=r;— 1 ' 



+C3 E \EnT)-Efm 

\a\=r]-l 
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Comparing the right hand sides in the estimates (3.4) and (3.7), we find 
similar terms involving the integral of 

(4.3) D\ = Ql^'^{TPu,T^u)dT'^u^{t), Df = hk{x)Ql^''{T^u,r^u)n^ (t) 

for |a| < m — 1, /3 + 7 + /X = Q and I7I < m — 1, 
(4.4) 

for \a\ = m — 1, 

(4.5) Dl = Q^K^^'dfidu^d^T'^u^d^T'^u^ , Z)| = 
for |a| = m — 1, and the semilinear terms 

(4.6) D\ = Sf^'^'iT^u, T^u)dT''u^{t), Dj = hk{x)Sj,'"{T'^u, T^u)T''u^{t) 

for |a| < ?Ti — 1, /3 + 7 + /X = a. 

To prove the claim, we need only to control these terms for m = k, 77 
separately. The proof will be similar to the proof of (8.5) and (8.6) in [46]. 
Here, we give the proof of the bound for the terms occurred on the right 
hand side of (3.7), that is, to control 



and 



hl{x)D]dx 



dUx 



with j = 1,2,3,4. The bound for the other terms follows by the same 
argument. In principle, the quasilinear and semilinear terms in Ij has been 
obtained in [46] and [8] separately. The new terms are the terms in //-,-. 
Here, we give the proof of both estimates of Ij and Ilj for completeness. 

4.1. Higher energy. For the first series of estimates we take m = k in 
(4.3)-(4.6). Recall that \hi\ < 1, \hkix)\ < r-V2(i + r)-i/2 ^nd using 
Hardy's inequality, we obtain immediately 



(4.7) 



< 



M 

C J2 Yj \\dTf^u^d'^T^u-^\\L2\\dT''u\\L2 

I,J=1 |o|<ft-l,/3+7<a,|7|<K-2 
M 

+^ Yl Yl \\dT^u^dT'^u-^\\L2\\dn\\L2. 

I,J=l |o|<K-l,/3+7<Q 

For the first term on the right-hand side of (4.7), we have either |/3| < k' or 
I7I < k'— 1, with k' = [k/2]. Note that since k > 9, we have k'+3 < k—2 = i]. 
We will also use that (t) < C{r){cjt — r). 
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In the first case, we estimate using (2.3) and (2.10b) 

\\dTf^u^d'^T"'u-^\\L2 <C{t)-^\\{r)dT'^u^\\L^\\{cjt - r)d'^T"'u'^\\L2 
<C(t)-ii?l^/l%(n(t))X,(n(t)) 
<C{tr'El/\u{t))El/\uit)). 

In the second case, we use (2.5) and then (2.10a) 

||9rV92rV||i2 <C{t)-^\\dT'^u\\L2\\{r){cjt - r)d^T^u'^\\L^ 
<C{t)-'El/\u{t))X^^^+,{u{t)) 
<C{ty'El/\u{t))X,{u{t)) 
<C{tr'El'\u{t))Ey\u{t)). 

For the second term on the right-hand side of (4.7), as in the proof of 
Lemma 2.6, we can use (2.6) and (2.5) to get (assuming |/3| < I7I) 

< C(t)-i(|Kc,t-r)arV(t)||ioo + ||(r)ar'5n^(t)||i^)||an^'^(t)||i2 

< C{tr\Ell\u{t)) + X,{u{t)))El/\u{t)) + CEy\u{t))El/\u{t)) 

< C{tr'Ell\u{t))Ell\u{t)) . 

Going back to (4.7), we have established the inequality (4.1). 

4.2. Lower energy. The second series of energy type estimates with m = rj 
will exploit the null condition. Let cq = min{c/}, the integrals will be 
subdivided into separate integrals over the regions r < cot/2 and r > cot/2. 
Inside the cones. On the region r < cot/2, we have that 

I,J,K |a|<»;-l,/3+7<a 



'U^ \\L^r< cot/2) ■ 



Since r < CQt/2, we have (t) < C{cit — r) for any I. Recall also + 3 < k, 
ItI + 2 < ?7, and |a| + 1 < ry, thus, using (2.4) and Lemma 2.7, the terms in 
the first part of the right hand side can be estimated by 

< C{ty'^''^\\{cit - r)'^l'^dT^u^\\L^\\{cjt - r)d'^T^u\\L2\\dT''u\\L2 

< C(t)-3/2[EV2^3(^(t)) + X|^|+3(n(t))]X|^l+2(n(t))^y2(^x(t)) 

< C{t)-^'^E,{u{t))El'\u{t)). 
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For the terms in the second part, using (2.6) instead of (2.4) and assuming 
|/3| < |7| (and so \/3\ < [t] - l]/2 < r/ - 3), we see 



< {t)-^/^\\r{cit-r)'/^dr''u'\\Lo 



{cjt 



< {tr'/' (i?|'/|%(n(t)) + X|^|+3(u(t))) (ei/^^m^)) 

< {t)~^/^E,{u{t))El/\u{t)), 



Ey\u{t)) 

,l+2(n(t))) Ey\u{t)) 



L2 

+ x 



where in the third inequality we have used the Hardy inequahty. 

The estimates for IIj{r < CQt/2) proceeds similarly with the obvious 
modification by using the Hardy inequality for terms involving r°n (and 
\hk{x)\ < C7r-V2(r)-V2 < c/r). This gives us the required upper bound 
for the portion of the integrals over r < CQt/2 in (4.2). 

Away from the origin. It remains to give the estimate for r > cot/2. 
It is here, finally, where the difference of speed cj and the null condition 
enters. 

Non-resonance. Let us start with non-resonant terms, that is, those for 
which (I, J,K) / {K,K,K). In this case, for Ij + IIj{r > cot/2), we need 
to control the quasilinear terms 
(4.8) 



{I,J)f^{K,K) [S+'y<a,\a\<ri-l,j<ri-2 

and the semilinear terms 
(4.9) 

E E 

{I,J)^{K,K) \a\<ri-l,l3+j<a 



'i/2(r)i/s 



L^r>cot/2) 



^1/2(^)1/5 



L^r>cot/2) 



We separate two cases for (4.8): I ^ J and I = J ^ K . In the first case, we 
have ci ^ cj, and 

(t)3/2 < C{r){{cit - r) + {cjt - r))i/2 < C{r){cit - r)^/'^{cjt - r)^'"^ . 

Using (2.4) and Hardy's inequality we have the estimate 

\\dTPu'd^T^u'dT^u''\\Li(,^,^,t/2) + \\dT^u'd^T^u-^T^u^^ 
< C(t)-3/2||(r)(c/t - r)^^^dTl^u^\\Lo^\\{cjt - ry/^d^T^u-^\\L2 
X (||5r°n^||i2 + ||r%^/r||^2) 



< C{t)'^/'[E^l^_^^{u{t)) + X\p\+3{u^^^^^ 

< C{tr^'^E,{u{t))El'\n{t)). 
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Otherwise, HI = J ^K, (2.4) and (2.2), we get 
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\\dT^u'd^T^u'dT''u^h^r>cot/2) 

+ ||r"V2^,)-l/25r/.^/52r7,/r%^||^.(,>,,,/2) 

< Ct'^/^t)'^\\dT^u^\\L2\\{cit - ry/^d^T^u^\\L2 

X (\\{r){cKt-r)^/^dr^u''\\L^ + |Kr)V2r°n^||^^) 

< Ct~'/'{ty^Ell^^^{u{t))X\^\+2Ht)) [El^^,{u{t)) + X|,|+3(n(t))' 

< Ct-y'{t)-'E,{u{t))El/\u{t)). 



For the semihnear terms (4.9), we have I ^ K and so c/ 7^ ck- Assuming 
ci < ck, then 



|r"u 



rV2(r)i/2 



L^{r>CQt/2) 



'u 



LI 



a/2r-^X| 



< Ct-V2(t)-ii?^(^(i))i^i/2(^(i)) 



by using (2.4) and (2.2). The case cj > ck can be handled the same way. 
This completes the estimates for non-resonant terms over r > cot/2 in (4.2). 

Resonance. In the resonant case, we have (/, J, K) = (K, K, K) and we 
will denote = u and ck = c. An application of Lemma 2.2 yields the 
following upper bound for Ij{r > cot/2): 



WTT^ud'^T^udT'' 



^+7<a,|7|<??-2,|al<»;-l 

+ \\dTf'udTT^udT''uhi^r>cot/2) 

+ 11 (ct - r)dT^ud^T^udT''u\\mr>cot/2) 

/3+7<a,|a|<f?— 1 

+ 11 (ct - r)ar'^n5nnar°u||ii(,>,„i/2) 



'l^\\L^(r>cot/2) 



rr^«an«5r-n^||ii(,>,„i/2) 



We still need to get an additional decay factor of t 
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Since r > cot/2, we have (r) > C{t). Thus, we have using (2.2) 

\\Tr^ud^T'yudT^u\\L^r>c,t/2) 

< C{t)-''^El'\u{t))E,{u{t)). 
In a similar fashion, the second term is handled using (2.3): 

\\dT^udTT'<udT^u\\L^t^r>c,t/2) 

< C{tr' \\dT^uh2 II {r)dTT^u\\L^(^r>cot/2) \\dT''uh2 

< C{tr'El/^^.-^{u{t))E,{u{t)) 

< C{tr'Ei'\u{t))E,{u{t)). 

The third term is estimated using (2.3) again and (2.10a). 

II (cf - r)dT^udT<udV''u\\L^t^,^,^t/2) 

< C(t)-i(r)5r^n||i^(,>,„,/2)|Kct-r)52nn||i2||9r"u|U2 

< C(t)-ii?|^/|'^3(u(t))X|^|+2(n(t))ii;y2(n(t)) 

< C{tr'Eil\u{t))E,{u{t)). 

For the terms arising in semilinear part, by (2.2), (2.6) and the fact (r) > 
C(t), we have 

||rr^nannar-n^||ii(,>,„i/2) + \\{ct - r)9r^n5nn9r°n||ii(,>,„i/2) 

< {ty^'^{W)^''^^^^AL'^ + \\{r)^''^ {ct - r)dT^u\\L^)\\dT^u\\L2\\dT''u\\L2 

< {tr^'\E\l^l^{u{t)) + X|^n.3)£;,(n(t)) 

< C{tr'I^Eil\u{t))E,{u{t)). 

To complete the proof of (4.2), we still need to give the estimate for IIj{t). 
In the resonant situation, noting that we always use norm to control the 
terms involving dT°^u in the proof of /j(t), it is easy to adapt the previous 
proof to get the required estimate for II j (t) by using Hardy's inequality and 
the fact that |/ifc(x)| < Cr~^. 



4.3. Conclusion of the proof. We are now ready to complete the proof 
of Theorem 1.1, by using (4.1) and (4.2). 
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Recall the definition (3.3) and (3.6) for E^, and E^j^, we know by 
Sobolev embedding, 

m{t)\+sMEl^{t)\) < C\\du\\L^EMt)) 

1 I 1 k>0 
\a\=m—l — 

< CEl/\u{t))EMt)) 

< C^CieEm{u{t)) 

for some C4 > 1, and 1 < m < k. 

Based on this observation and the smallness assumption of e and 6 (such 
that 8C3CiCi{e + 5) < 1), we can easily obtain the following inequalities 
from (4.1) and (4.2), 

(4.10) LE^{T)+ sup E^{u{t)) 

0<t<T 



< 2C3^,(n(0)) + 2C3 / {t)-'El/\u{t))E^{u{t))dt 

Jo 



and 



(4.11) LE^{T)+ sup E^{u{t)) 

0<t<T 



< 2C3E^{u{0)) + 2C3 r^'^{t)'^El/\u{t))E^{u{t))dt 

Jo 

on any interval [0, T] with 

T<To = sup{r : ^y2(n(t)) < 2Cie,t G [0,r]} . 

1 /2 

Since Erj {u{t)) < 2Cie, an application of the Gronwall inequality to 
(4.10) gives us 

LE^it) + E^{u{t)) < 2C3e4^^^i^ii;,(n(0))(t)^^^^i^ 

Inserting this bound into (4.11), we obtain 

LEr,{t) + Er,{u{t)) < 2C3Er,iu{0)) 

X exp (^(2C3)3/2e2^3C,.^l/2(^(o))^"(t)-^+2C3Cie^-l/2^^ 

Setting Ci = ^/2C^, 



and C2 = max(Cf e^*^!^, 2Cf ), then we have 

LEr^it) + E^{u{t)) < 2C3^^(n(0))exp(2Co^y2(^(o))) < 2C3e2 = [Cisf , 

LE^it) + E,{u{t)) < C2E^iu{0)){tf^' . 

To ensure the finiteness of Co, in addition to the smallness assumption on e 
to absorb the perturbation, we need also require Cfe < 1/2. 
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1 /2 

With this we see Erj {u{t)) remains less than Cie throughout the interval 
< t < Tq. a standard continuity argument shows that Eri{u{t)) is bounded 
for all time, which completes the proof of Theorem 1.1. 

5. Appendix: the case of asymptotically flat manifolds 

In this appendix, we will give the almost global existence (and global 
existence respectively) for the Cauchy problem of the quasilinear wave sys- 
tem (1.3) when spatial dimension is three (and higher), posed on certain 
asymptotically flat manifolds. 

Let us begin with the space-time manifolds. We consider the asymptoti- 
cally flat Lorentzian manifolds (M^'''"',^) with 

n 

Q = Qaisit, x) dx°' dx^ = Qapit, x) dx"' dx^ . 

Q,/3=0 

The metric g is assumed to be a small asymptotically flat perturbation of 
the Minkowski metric. More precisely, we suppose gai3(t,x) € C°°(M'^"^") 
and, for some fixed p > and 5 <^ 1, 

(H2) V7 G Ni+" \dl,{g^^{t,x) - m„^)| < 

with {niais) = Diag{ — 1, 1, 1, • • • ,1) being the standard Minkowski metric 
and (x) = Y^l + l^p. An example of such metric can be 

9ap = rua^ + S{xyP + 5^{t/{x)){x)-P 

with (f) G Cq°. Since (5 <C 1, it is clear that the metric g is a non-trapping 
perturbation. Let g = — det{gai3), the Laplace-Beltrami operator associated 
with g is given by 

□g = Vg-'dag'^^^d/,, 

where {g°'^{t,x)) denotes the inverse matrix of {gaf3{t,x)). 

We would also like to investigate the case of radial metric, by which we 
mean that, when writing out the metric in polar coordinates {t,x) = {t,roj) 
with oj G we have 

5 = gooit,r)dt^ + 2goi(t,r)dtdr + gii{t,r)dr'^ + g22{t,r)r'^duj^ . 

In this form, the assumption (112) on asymptotically flatness is equivalent 
to the following requirement 

(H2') \dl,{~goo + 1,511 - 1,522 - l,5oi)| < C^S{x)-^^^-P. 

Consider the initial value problem for the quasilinear wave equations of 
the form 

(5.1) {nguY = N^{u,u),I = 1,2,- ■■ ,M 

in which the quadratic nonlinearity iV = Q + S is of the form (1.4). Our 
construction of solutions will depend on the energy integral method, which 
requires the quasilinear part to be symmetric (1-5). 
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In contrast to the null- form system, as in [41], we will be able to avoid 
the use of the scaling vector field S, and the vector fields to be used will be 
labeled as 

Y = {Yo,...,Y^) = {d,n). 
For the energy norm, we will use the standard energy norm 

1 ^ r 

Ei{uit)) = -Y, / \du'{t,x)fdx. 
1=1 -^^^ 

The higher order derivatives will be estimated through 

(5.2) Em{u{t))= E EiiY'^uit)), m = 2,3,--- . 

\a\<m—l 

As for the null-form systems, an important intermediate role will be played 
by the following local energy norm 



AI 



1=1 |a|<m,-l 



r"u-f(i)| 



with fj, G (0, 1/2) and fi' > fi to be determined later. This norm is extracted 
from the local energy norm, which is defined as 

ft 



(5.3) 



LEm{t) 



'Jmiu{T))dT . 



In the case of fi' = fi, such norm will be the KSS norm and we will denote 
such norm by KSSm.{i)- We will choose fi = 1/4. When n > 4, the choice 
for n' will be min(n — 2, 2p — 1, 3)/4 for the general case and n' = min(n — 
2, 1 -|- 2/3, 3)/4 for the radial metric. In the case of n = 3, we will set 
fi' = min(2/? — 1, 3)/4 (and n' = min(2p -|- 1, 3)/4 for the radial metric). 
In order to describe the solution space, we introduce 

i/f (M") = {/ G L\W) : y"/ G L^, \a\ < m}, 



with the norm 
(5.4) 



\a\<m 



Solutions will be constructed in the space HyiT) obtained by closing the 
set C7°°([0,r);C7^(M")) in the norm sup EU^{u{t)). Thus, 

0<t<T 



m—l 



L{t,x):du{t,-)e fl C^([0,r);/7™-i-^) 
j=0 



By Sobolev embedding, it follows that ijp(r) C C™-["+2]/2qq^ ^ M"). 
Let us now state our main result precisely. 
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Theorem 5.1. Let n>3, 6<^1, p>l for the general metric and p > 
for the radial metric. Assume that the nonlinear terms in (5.1) satisfy the 
symmetric condition (1.5). Then there exist constants eQ,co <C I, such that 
the Cauchy problem for (5.1) has a unique global (almost global for n = 3) 
solution u e Hyi^) /^^^ * ^ ^^'^ 

^ ^ f oo n > 4 

^ I exp(co/e) n = 3 ' 

when the initial data satisfy 

(5.5) ^y2(u(0)) = e < eo,K > n + 4. 

Moreover, the solution satisfies the bounds for some Ci > 1, 

sup E^{u{t)) + LE^{Te) + 5-ini^KSS,{Te) < Cy. 

t(^[0,T,] 

5.1. Commutation with vector fields. In preparation for the energy es- 
timates, we need to consider the commutation properties of the vector fields 
Y with respect to the nonlinear terms. 

Lemma 5.1. Let u be solution of (5.1). Assume that the nonlinearity is of 
the form (1.4). Then for any a G N'', 

/3+7+At=o |/9|<|a|-l 

in which each is a quadratic nonlinearity of the form (1.4), and r^ with 
m € N denote functions such that 

|9°r„(t,x)| < C„5(r)-^-"^-l°l for any a G N^+" . 

Moreover, if |/i| = 0, then N'^ = N . In addition, if the metric is radial, 

□gy"u= Y N;;[{y^u,Y^u)+ Y {riV^T^u + r2VT^uy 

Proof. It is easy to check that 

[Y, N] {u, v) = YN{u, v) - N{Yu, v) - N{u, Yv) 

is a quadratic nonlinearity of the form (1.4). 
By (H2), we have 

□0 = + rod"^ + rid , 
with □ = Diag{-d'^ + A, • • • , -d^ + A). Recall that 

[□,y,] = o, 

we want to prove the result by induction. It is clear the result is true for 
\a\ =0. Now assume that it is true for any a with |q;| = m. Given with 
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jaol = m + 1, we can find some j and a with \a\ = m and y"'' = YjY". 
Then by the inductional assumption, we can calculate as follows 

□gy°o^ = C\gYjY°'u 

+ YjNl^{Y^u,Y"'u)+ Yj (^rod'^Y^u + ndY'^u 

= Y (rod^Y'^u + ridYf^u^ 

+ Yl {[Yj,N;^]{Y^u,Y^u) + N;^{YjY'^u,Y^u) + N;^{Y'^u,YjY^u)^ 

/3+7+/i=a 

which is of the required form. This completes the proof for the general case. 
When the metric is radial, the rotational vector fields Q are commutative 
with Dg and we need only to give the estimate for d. In this case, we have 
[r2~jd^,d] = r^-jd^ with j = 1,2 and the same argument will give the 
proof. □ 



5.2. Sobolev-type inequalities. As in Metcalfe-Sogge [41], we need to use 
the Sobolev inequalities which does not involve the Lorentz boost operators 
and scaling vector field, which is also related with the trace estimates. 

Recall that we have the following trace estimate (see (1.3) of Fang and 
Wang [7]) 

r^||/(ro.)|| < C||/||^„/,_.,0 <s<{n- l)/2 

which gives us (with s = 1/4) 

(5.6) r'/'\\f{ru:)U^<C ^ Wd^n^fh^- 

\a\ + \l3\<{n+3)/2 

Moreover, Sobolev embedding in the polar coordinates gives us for r > 2 
(see e.g. (2.1) in Metcalfe-Sogge [41]) 

(5.7) A^-'^^WfirOjU^ <C Y l|5"^^/llL2(r-l<|x.|<r+l) • 

la| + l/3|<n/2+l 



28 CHENGBO WANG AND XIN YU 

As a direct consequence of the estimates (5.6) and (5.7) is 

(5.8) WfgWh 

< \\f9\\ 

i>3 

< llr-V4f||2 ||r-l/4„||2 

^ If /IIl°°(|x|<2)IF fi'llL2(|2,|<2) 

+ V llr("-l)/2i.||2 , , , Jlr-("-l)/2 l|2 

+ lli°°(i-l<|a;|<i+l)IF 3\\L2{j-l<\x\<j+l) 

la|<(n+3)/2 

+CY1 Yl Il^".^lli2(i-2<|x|<j+2)lk ^" ^^^^fi'llL2(j-l<|x|<i+l) 
i>3 |Q|<(n+3)/2 

< C||r-l/4(r)-("-2)/45||2, ^ ||r-V4(^)-{-2)/4ya^||2^ 

|a|<(n+3)/2 



5.3. Energy and local energy estimates. As in the beginning of Section 
3, assume that u{t) G HyiT) is a local solution of the initial value problem 
for (5.1). Our task will be to show that Ei^{u{t)) remains finite for alH > 
when n > 4 (and for some interval [0, exp(c/e)] with c <C 1 when n = 3). To 
do so, we will derive integral inequalities for E^{u{t))+LEf^{t). If (5.5) holds, 
then Sy^(u(0)) = e and E^{u{t)) + LE^{T) + (log(2 + T))-^KSS^{T) < 
(Cie)2 for certain small interval t G [0,T]. Define 



To = sup{r : E^{u{t))+LE^{T)+{log{2+T)r'KSS^{T) < {2Cie)\t £ [0,r]} . 



Here the constant Ci > 1 will be determined later. All of the following 
computations will be valid on the time interval [0, Tq). 

It is easy to check that the same argument of Section 3 gives us the 
required local energy estimates, that is, by writting 



□g = □ + (5°^ - m"^)a„5;3 + nd , 
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and applying Lemma 2.4 with /i"'^ = g"'^ — m"^, Lemma 5.1 and the sym- 
metry condition (1.5), we have 

sup EmMt)) + LEm(T) + (log(2 + T))-'KSSm.{T) 

< CEM0)) + C6LEm.{T) 



|«| <m— 1 



T 



C^{t)dt 



+ c sup 

^ fe>0 



C?,it)dt 



\a\<m~l 

< CEM^)) + CCi{5 + e)LEm{T) 

+c J2 \EnT)-Eno)\+csup iei^t) - E^^.m 

\a\=m—l 
M 



A:>0 



a\=m — l 



sup 

fc>0 



1=1 |a|<m — l,/3+7+/i=a 
M 

+CE E 

1=1 \a\<m—l,l3+'y+fi=a,\'y\<m—l 
M 



sup 

fc>0 



/ / Sl'''(rl^u,T^u){dt,Lk)T''u^{t)dxdt 
Jo Jr^ 

r [ Ql^''ir^u,r^u)idt,Lk)r''u'it)dxdt 

Jo Jr^ 

+'^fl Yl ^^^P / Qji^^d^dpu'^df,T''u^{dt,Lk)T''u^dxdt 

I I , fc>0 Jo 
-' = -1- |a|=m— 1 ~ 

M T I- 

+'^1^ / / QjK'^{dt,hidi)d,3u-^df,r'^u^d^r'^u^dxdt , 

I I , fc>0 JO Jr^ 
J = l |a|=m— 1 — 

where we have used the assumption 

(5.9) 2^' < p - 1 + 2^, p > 1 

for general metric. When the metric is radial, the assumption can be weak- 
ened to be 



(5.10) 



2/u' < p + 2fi, p> {g radial) . 



Here, as before, Lk = fk{dr + ]:) = hl{x)di + hk{x) with hk{x) = fk{\x\) /\x\ 
and /i*^(x) = x*/ifc(x), /i = 1/4 < /i' < 3/4, 

M 

(5.11) E^{t) = y^iQ'j^d^, - -QYD / d0'd^T^u''d,T^u'dx, 

1=1 J^^ 

and 

M . 

(5.12) £;2^ = VQj'J'^" / dpu'd^T^u''Lkn'dx. 
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5.4. Conclusion of the proof. We are now ready to complete the proof of 
Theorem 5.1, by using (5.8) and the local energy estimates in the Subsection 
5.3. 

We need only to give a better upper bound for 

sup E^{u{t)) + LE,{T) + (log(2 + T))-^KSS^{T) . 

By the local energy estimates in the Subsection 5.3, we see that there exists 
an universal constant C2 > 1 such that 

sup E^{u{t)) + LE^{T) + (log(2 + T)r^KSS^{T) 
< C2E^{um + C2Ci{5 + e)LE^{T) + C2 ^ ^MEl^T) - E^^^ 



|a|=K— 1 



fc>0 



+C2 Yl \E?{T)-Ef 



|o|=K — 1 



|q|<k-1 \ 



r 



|al<K-l,l/3|<(K-l)/2 



dy°'u\\L^Ll{lO,T]xR") + 



< C2E^{um + C2C^{5 + e)LE^{T) + C2 ^M^lk^T) - E^^m 



\a\=K-l 



k>0 



+C2 Y \EnT)-Eno)\ 

\a\=K—l 

+ C2 sup V ||^-l/4^^^-(n-2)/4gya^ 

te[o,T] 



Il2l2([0,T]xR")' 



\a\<K-l,{K+n+2)/2 

where in the last inequality, we have used Hardy inequality and (5.8). 

Recall the definition (5.11) and (5.12) for Ef and E^j^, we know by 
Sobolev embedding, 

Y {\Ent)\+sup\El,m < CWduU^EMt)) 

\a\=m~l — 

< CEl'\u{t))EMt)) 

< CsCisEMt)) 

for some C3 > 2. 

Based on this observation and the smallness assumption of e and 6 (such 
that 8C2C3Ci(e + S) < 1), we can easily obtain the following inequality, 
when (k + n + 2)/2 < K — 1, i.e. k > n + 4, 

(5.13) sup E,iu{t)) + LE^{T) + (log(2 + T)y'KSS,{T) 

ie[o,r] 

< 2C2^«(u(0)) +2C2i5;y'(log(2 + T)f'-{LE^{T) + (log(2 + T))-^ K S S ^{T)) 
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on any interval [0, T] with 



T < To = sup{r : E^{u{t))+LE^{T)+{log{2+T))-'KSS^{T) < {2Cief,t€ [0,T]} . 



(5.14) sup E^{u{t)) + LE^{T) + (log(2 + T))-'KSS^{T) < SCfe^ 



by setting Ci = \/^2 and 5, e < eo = l/(16C3C2Ci). In the case of n = 3, 
we have (5.14) for T < and T < Tq, if we set Ci = \fC^^ <5, e < eo = 
l/(16C3C2Ci), and 



with Co = l/(32C2Ci). When n > 4, we denote = oo. Here, we may need 
to let b even smaller such that the assumption in Lemma 2.4 is satisfied. 

1/2 

With this we see the sum of E^ {u{t)) and local energy norms remains less 
than \/3Cie throughout the interval < t < To (and t < T^). A standard 
continuity argument shows that Ei^{u{t)) is bounded for t S [0,Ts), which 
completes the proof of Theorem 5.1. 
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